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Abstract. Modifications of photoionization arising from resonant electron-electron 
correlations between neighbouring atoms in an atomic sample are studied. The 
Y^ . sample contains atomic species A and B, with the ionization potential of A being 

•'-H I smaller than the energy of a dipole-allowed transition in B. The atoms are subject 

^»L/ to an external radiation field which is near-resonant with the dipole transition in 

r^ , B. Photoionization of an atom A may thus proceed via a two-step mechanism: 

photoexcitation in the subsystem of species B, followed by interatomic Coulombic 
decay. As a basic atomic configuration, we investigate resonant photoionization in a 
three-atomic system A-B-B, consisting of an atom A and two neighbouring atoms 
B. It is found that, under suitable conditions, the influence of the neighbouring 

fT^ . atoms can strongly affect the photoionization process, including its total probabilty, 

^^ \ time development and photoelectron spectra. In particular, by comparing our results 

with those for photoionization of an isolated atom A and a two- atomic system A-B, 

v_j , respectively, we reveal the characteristic impact exerted by the third atom. 
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^- 1. Introduction 



Photoionization of atoms and molecules is one of the most fundamental quantum 
processes. It played a key role in the early days of quantum mechanics and has ever since 
been paving the way towards an improved understanding of the structure and dynamics 
of matter on a microscopic scale. Today, kinematically complete photoionization 
experiments allow for accurate tests of the most sophisticated ab-initio calculations. 
Besides, photoionization studies in a new frequency domain are currently becoming 
feasible by the availability of novel xuv and x-ray radiation sources [U l2l E] , giving rise 
to corresponding theoretical developments (see, e.g., [HEllG]). 

Various photoionization mechanisms rely crucially on electron-electron correlations. 
Prominent examples are single-photon double ionization as well as resonant 
photoionization. The latter proceeds through resonant photoexcitation of an 
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autoionizing state with subsequent Auger decay. In recent years, a similar kind of 
ionization process has been studied in systems consisting of two (or more) atoms. Here, 
a resonantly excited atom transfers its excitation energy radiationlessly via interatomic 
electron-electron correlations to a neighbouring atom leading to its ionization. This 
Auger-like decay involving two atomic centers is commonly known as interatomic 
Coulombic decay (ICD) [3 [8]. It has been observed, for instance, in noble gas dimers 
and water molecules [9]. In metal oxides, the closely related process of multi-atom 
resonant photoemission (MARPE) was also observed |10j . 

We have recently studied resonant two-center photoionization in heteroatomic 
systems and shown that this ionization channel can be remarkably strong [HI [121 US] ■ 
In particular, it can dominate over the usual single-center photoionization by orders 
of magnitude. Besides, characteristic effects resulting from a strong coupling of the 
ground and autoionizing states by a relatively intense photon field were identified. Also 
resonant two-photon ionization in a system of two identical atoms was investigated [H] . 
We note that photoionization in two-atomic systems was also studied in [T5l [T6] and 
pTllTS]. The inverse of two-center photoionization (in weak external fields) is two-center 
dielectronic recombination |19]. 

In the present contribution, we extend our investigations of electron correlation- 
driven interatomic processes by considering photoionization of an atom A in the presence 
of two neighbouring atoms B (see figure 1). All atoms are assumed to interact with each 
other and with an external radiation field. We show that the photoionization of atom 
A via photoexcitation of the system of two neighbouring atoms B and subsequent ICD 
can be by far the dominant ionization channel. Moreover, we reveal the characteristic 
properties of the process with regard to its temporal dependence and photoelectron 
spectra. In particular, by comparing our results with those for photoionization in a 
system of two atoms A and -B, we demonstrate the influence which the presence of the 
second atom B may have. 
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Figure 1. Schematic illustration of photoionization of an atom A in the presence of 
an external laser field and two neighbouring atoms B and B' . Apart from the direct 
photoionization of A there are interatomic channels via resonant photoexcitation of 
the "molecular" system B-B' and subsequent ICD. 
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Atomic units (a.u.) are used throughout unless otherwise stated. 

2. Theoretical Framework 

Let us consider a system consisting of three atoms, A, B and B' , where B and B' are 
atoms of the same element and A is different. We shall assume that all these atoms 
are separated by sufficiently large distances such that free atomic states represent a 
reasonable initial basis set to start with. 

Let the ionization potential I a of atom A be smaller than the excitation energy ^Eb 
of a dipole- allowed transition in atoms B and B' . Under such conditions, if our system is 
irradiated by an electromagnetic field with frequency uq ~ AEb, the ionization process 
of this system (i.e., essentially of the atom A) can be qualitatively different compared 
to the case when a single, isolated atom A is ionized. Indeed, in such a case A can 
be ionized not only directly but also via resonant photoexcitation of the subsystem of 
B and B', with its consequent deexcitation through energy transfer to A resulting in 
ionization of the latter. 

In the following, we consider photoionization in the system of atoms A, B and 
B' in more detail. For simplicity, we suppose that the nuclei of all atoms are at rest 
during photoionization. Denoting the origin of our coordinate system by O, we assume 
that the nuclei of the atoms B and B' are located on the Z-axis: R^ = (0, 0, Zb) 
and Hb' = (0,0, Z^/). The coordinates of the nucleus of the atom A are given by 
R^ = {Xa, Ya, Za). The coordinates of the (active) electron of atom A with respect to 
its nucleus are denoted by rA, where A G {A, B, B'}. 

The total Hamiltonian describing the three atoms embedded in an external 
electromagnetic field reads 

H = Ho + Vab + Vab' + Vbb' + Wa + Wb + Wb^, (1) 

where Hq is the sum of the Hamiltonians for the noninteracting atoms A, B and B'. 

We shall assume that the (typical) distances AR between the atoms are not too 
large, AR ^ c/AEb, where c is the speed of light, such that retardation effects in the 
electromagnetic interactions can be ignored. If transitions of electrons between bound 
states in atoms B and B' are of dipole character, then the interaction between each pair 
of atoms (A, 7) (with A, 7 G {A, B, B'}) can be written as 

where Ra,7 = Ra ~ R7 and 6ij is the Kronecker symbol. Note that in ([2]) a summation 
over the repeated indices i and j is implied. 

In ([2]), W\ denotes the interaction of the atom A with the laser electromagnetic 
field. The latter will be treated as a classical, linearly polarized field, described by the 
vector potential A(r, t) = Aocos(a;ot), where Aq = cFq/uq, Uq = cko is the angular 
frequency and Fq is the field strength. The interaction Wx then reads 

Wx = ^A{rx,t)-px, (3) 
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where pa is the momentum operator for the electron in atom A. 

Our treatment of photoionization will be based on the following points: 
Oscillator strengths for dipole-allowed bound-bound transitions can be very strong. 
This means that, provided that the distances between all the atoms in our system are 
of the same order of magnitude, the interaction between atoms B and B' is much more 
effective than the interaction between atoms A and B (or A and B'). Besides, atoms B 
and B' will, in general, couple much more strongly to a resonant laser field than atom 
A. In what follows, we shall assume that the intensity of the laser field is relatively 
low such that the interaction between atoms B and B' changes the states of the system 
more substantially than the coupling of these atoms to the laser field. Therefore, we 
shall begin with building states of the B-B' subsystem in the absence of the field. The 
second step of our treatment will be to include the interaction of the B-B' subsystem 
with the laser field and, in the third step, we complete the treatment of ionization by 
considering the interaction of atom A with both the laser field and the field-dressed 
subsystem of atoms B and B'. 



I. We denote the ground and excited states of the undistorted atoms B and B' by 
(po, 0e and 0Q, 0g, respectively. Let the corresponding energies of these states be Eq and 
Ee- The state ipBB' of the B-B' subsystem can be expanded into the "complete" set 
of undistorted atomic states represented by the configurations (i) 0o0o; (ii) 0o0e5 (iii) 
0e0o and (iv) 0e0e- In the approximation, which neglects the interatomic interaction, 
the configurations <;/>o0e ^^^ 0e0o are characterized by exactly the same value of the 
(undistorted) energy Eoe = ^o + ^e- The latter, in turn, strongly differs from the 
energies Eqq = 2eq and E^e = 2ee which are characteristic for the configurations 0o0o 
and 0e0e; respectlvcly. Therefore, provided that the distance between the atoms is not 
too small, the interaction Vbb' "will strongly mix the configurations (ii) and (iii) only, 
while the other configurations (i) and (iv) will be affected only very weakly. Taking this 
into account, it is not difficult to find the states of the subsystem of interacting atoms 
B and B' which read 

Vg = 0000 

V^+ = ^ (0e0O + 0O0e) 



1 



' %€) 



¥'- = -71 [(Pe(Po 

Ve = M'e- (4) 

These two-atomic states are normalized and mutually orthogonal. They posses energies 
given by Eg = 2eo, -E+ = Eo + Ee + vbb', E_ = E^ + Ee — vbb' and Ee = 2Ee, respectively, 
where vbb' = \4>e4>'o ^bb' 4'o4>'e)- Note that, for definiteness, vbb' has been assumed 
to be real and negative here, as will always be the case in our examples below (see 
section 3). 
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II. Let us now consider two interacting atoms B and B' embedded in a resonant 
laser field. One can look for a state of such a system by expanding it into the new set 
of states given by Eq. (II]), 

i){t) = g{t)ipg + a+{t)ip+ + a^{t)ip_ + 6(t)v?e- (5) 

Inserting the expansion (|5]) into the corresponding wave equation, we obtain a set of 
coupled equations for the unknown time-dependent coefficients g{t), a+(t), a_(t) and 
bit): 
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The system of equations iQ can be greatly simplified by noting the following. First, 
all transition matrix elements of the interaction with the laser field, which involve the 
asymmetric state </)_, are equal to zero and, thus, only the remaining three states can 
be coupled by the field. Second, if we suppose that the frequency of the laser field is 
resonant to the transitions ipg < — > ip+ and that the field is relatively weak such that 
the non-resonant transitions ip+ < — > ipe are much less effective than the above resonant 
ones, the system (E]) effectively reduces to 
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which can be readily solved by using the rotating wave approximation. Assuming that 
the field is switched on suddenly at t = 0, we obtain two solutions 
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In the above equations, we have introduced 

zi = - {Eg + E+-UJO- VLr) 

Z2 = ]^ {Eg + E+-UO + Qr) , (10) 



where Q,r = y{E+ — Eg — uJoy + 4\Wg^+\ is the Rabi frequency, Wg^+ — 
{ifg |Fo ■ {VB + PB') /(2a;o)| </'+) and W+,g = (W,,+)*. 

The two solutions in ^ correspond to two different initial conditions: at t = the 
system is either in the state (pg or in </?+. They are orthogonal to each other and form 
a "complete" set of field-dressed states of the subsystem B-B' . Note also that we have 
neglected the spontaneous radiative decay of the excited state </?+ which, in our case, is 
justified as long as |M/g^+| ^ F^, where V^ is the radiative width of v?+- 



III. Now, as the last step, we shall add atom A to our consideration. Let Xo and Xp, 
where p is the electron momentum, be the ground and a continuum state of a single, 
isolated atom A. The wavefunction of the total system A-B-B' can be expanded into 
the following "complete" set of states 

^{t) = ao{t) ^1 xo + /3o(t) ^2 Xo + y o?^pap(t) ^iXp + j d^pPp{t) ^2 Xp- 

(11) 

Here, the initial conditions are given by q;o(0) = 1, /3o(0) = and ap(0) = /3p(0) = 0. 
The coupling of atom A to both the subsystem B-B' and the laser field involves bound- 
continuum transitions which are normally much less effective than the bound-bound 
ones. For this reason, we may assume that the interactions of A with the laser field and 
the S-i?'-subsystem is weak and consider ionization of atom A in the lowest order of 
perturbation theory in these two interactions. As a result, by inserting the expansion 
flTTj) into the corresponding Schrodinger equation we obtain 

i-TT - ep«p = exp{-ieft) (i>i Xp \Wa + Vab + Vab'| ^1 Xo) 

'^ - </5p = exp(-ze^t) (v-a Xp \Vab + Vab\^i Xo) , (12) 

where e^ is the energy of the electron in the initial state xo of atom A and e^ is the 
electron energy after the emission. The probability for ionization of the three-atomic 
system, as a function of time, then reads 



P{t)= J d'p(\a^{t)\' + \^^{t)\'). (13) 

Note that equations ( !T2|) are readily solved analytically. However, the resulting 
expressions are somewhat lengthy and will not be given here. 
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Figure 2. Photoionization probability for Li, Li-He and Li-He-He systems in 

an external electromagnetic field, given as a function of time. The field strength is 
Fq = 10~^ a.u., the field is linearly polarized and its frequency is resonant to the 
corresponding transition in the He or He-He subsystem. The distances between Li and 
each of the He atoms is always 14 a.u. The atomic positions are aligned along the field 
polarization with the Li atom in the middle of the three-atomic system. The solid, 
dash and dot curves display results for Li-He-He, Li-He and Li systems, respectively. 
Note that the ionization probability for an isolated Li atom has been multiplied by a 
factor of 500. For more explanation see the text. 



3. Results and Discussion 



Based on the results obtained in the previous section, let us now turn to the discussion 
of some aspects of photoionization in a system consisting of one lithium and two 
helium atoms. We suppose that in our three-atomic system the positions of the 
lithium and helium atoms are given by the vectors Rlj = (0,0,0), Rhc = (0,0, Z) 
and Rhc' = (0, 0, —Z), respectively. Our system is initially (at time t = 0) in its ground 
configuration and is irradiated by a monochromatic laser field. The field is linearly 
polarized along the Z-axis and its frequency is resonant to the (pg - </?+ transition in the 
He-He subsystem, i.e., £"+ — Eg — uq = Q. 

Choosing Z = 14 a.u. we obtain that the energy spitting AE± = \E^ — E_\ 
between the states </?+ and if- of the He-He subsystem is 5.4 x 10~^ eV. Assuming 
a field strength of Eq = 10~^ a.u., the corresponding Rabi frequency amounts to 



Q 



^Hc 



2\W, 



1.3 X 10 ^ eV which is much less than AE±. 



In figure [21 we present the probability for ionization of our system as a function 
of time. The probability shows a non-monotonous behaviour in which time intervals. 
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when the ionization probabihty rapidly increases, are separated by intervals, when 
the probability remains practically constant, reflecting oscillations of the electron 
populations with the Rabi frequency ^r^^_^^ between the ground and excited states 
of the He- He subsystem in a resonant electromagnetic field. 

For comparison, we also show in figure |2]results for ionization of a single (separated) 
Li atom and for ionization in a two-atomic Li-He system. In the latter case, the lithium 
atom is located at the origin (Rli = (0,0,0)) and the coordinates of the helium atom 
are Rne = (0,0, 14 a. u.). The frequency of the laser field is assumed to be resonant to 
the ls'^^S-ls2p^P transition frequency of the corresponding bound states of a single 
He atom. 

In contrast to the single-atom ionization, in both the two- and three-atomic cases 
the ionization probability demonstrates a step-wise temporal development in which 
time intervals of rapid probability growth are followed by intervals of almost constant 
probability. We point out that in the three-atomic case, however, the size of these time 
intervals is shorter by a factor of a/2. 

Compared to ionization of a single Li atom, ionization in the two-atomic system 
is very strongly enhanced [HI [121 [E]- When the three-atomic system is irradiated, 
the enhancement increases even further. In the range of small values of t, where all 
ionization probabilities still increase monotonously, this additional enhancement is equal 
to a factor of 4. At larger t, however, when the two ionization probabilities exhibit step- 
wise behaviours, this additional enhancement due to the presence of the second He atom 
is reduced to a factor close to 2 on average, as can also be seen in figure O 

All the above features can be understood by noting the following: 

i) For the chosen set of parameters of our two- and three-center systems, the indirect 
channels of ionization, which involve two- or three-atomic correlations, are substantially 
stronger than the direct one. Therefore, these correlations have a dominating effect on 
the ionization. 

ii) At small t, ionization in the two- and three-atomic systems is basically a two-step 
process: the first step is photoexcitation in the He or He-He subsystem and the second 
step is a consequent energy transfer to Li. In each case, both these steps are described by 
basically the same dipole transition matrix element of the subsystem. Since, compared 
to a single He atom, this dipole element in He- He is by \/2 larger than in He, one obtains 
a factor of 2 for the enhancement in the ionization amplitude, leading to a factor of 4 
in the ionization probability (see also [20]). 

iii) At larger t, when Rabi oscillations show up, the second step "saturates" in the 
sense that the averaged probability to find the corresponding subsystem in the excited 
state becomes equal to 50%. Therefore, the ionization probability in the three-atomic 
system is now larger (on average) by a factor of 2 only. 

iv) The origin of the step-wise behaviours of the ionization probabilities for the 
two- and three- atomic systems lies in the oscillations of the population between the 
ground and excited states in the He atom (for the two-atomic case) or in the He-He 
subsystem (for the three-atomic case). The scale of these oscillation is set by the Rabi 
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Figure 3. Energy spectra of the emitted electrons, as a function of A = e^ — Cg — wq, 
for the same parameters as in figure^] The pulse duration is 100 ps. a) Solid and dash 
curves show results for ionization of Li- He- He and Li- He systems, respectively, b) Solid 
and dot curves display results for ionization of Li-He-He and Li systems, respectively. 
The results for the Li system have been multiplied by a factor of 500. 



frequency and, because in the He-He subsystem the latter is larger by a factor of a/2, 
the corresponding time intervals are shorter by the same factor. 

Additional information about the ionization process can be obtained by considering 
the energy spectrum of emitted electrons. Such a spectrum is shown in figure [3] for the 
same systems and parameters as in figure [2] and for a pulse duration of T = 100 ps. 

In panel (a), we compare the energy spectra of electrons emitted in the process 
of photoionization of Li-He-He and Li-He systems. In both cases, the main feature 
is the presence of three pronounced maxima. The origin of these peaks is similar to 
the splitting into three lines of the energy spectrum of photons emitted during atomic 
fluorescence in a resonant electromagnetic field [21]. In such a field, the ground and 
excited levels of the He and He- He subsystems split into two sub-levels, which differ by 
the corresponding Rabi frequency Vt^. As a result, the resonant electronic correlations 
between these subsystems and the Li atom lead to an energy transfer to the Li which 
peaks at ujq and ujq ± ^r/2. Since, as was already mentioned, the Rabi frequencies 
of these subsystems differ by a factor of \/2, the magnitude of the separation between 
the corresponding maxima in panel (a) of figure [3] also differs by this factor. Note 
also that the widths of these main maxima as well as the appearance of additional 
multiple maxima, seen in the figure, are related to the finiteness of the pulse duration; 
the distance between the latter is roughly given by 2tt/T. 

The distinct influence, which the interatomic electron-electron correlations exert on 
the shape of the photoelectron spectra, is further highlighted in panel (b) of figure [31 It 
compares the energy spectra of photoelectrons emitted from our Li-He-He system and 
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an isolated Li atom. In the latter case, there is only one main maximum, while the two 
main side peaks are missing, as one would expect (the additional multiple maxima are 
related again to the finiteness of the pulse duration). 

4. Conclusion 

We have studied resonant photoionization in a system A-B-B' consisting of three atoms, 
with two atoms B of the same element and one different atom A. We have shown that 
the mutual correlations among the atoms can largely enhance the ionization probability 
and distinctly modify also other properties of the process in a characteristic manner. In 
particular, as compared to the case of resonant photoionization in a two-atom system 
A-B, it has been demonstrated that the presence of a second atom B can (i) further 
enhance the photoionization process, (ii) change the time dependence of the ionization 
probability and (iii) move the side peaks in the photoelectron spectrum further apart. 
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